A Hybrid Density Functional Study of Oligothiophene/ZnO Interface for Photovoltaics 
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Organic/inorganic donor- acceptor interfaces are gaining growing attention in organic photovoltaic 
applications as each component of the interface offers unique attributes. Here we use hybrid-density 
functional theory to examine the electronic structure of sexithiophene/ZnO interfaces. We find that 
interfacial molecular orientations strongly influence the adsorption energy, the energy level align- 
ment, and the open circuit voltage. We attribute the orientation dependence to the varied strength 
of electronic coupling between the molecule and the substrate. Our study suggests that photovoltaic 
performance can be optimized by controlling the interfacial design of molecular orientations. 



Organic photovoltaics (OPV) differ from conventional 
solar cells in that their fundamental mechanism is dom- 
inated by interfacial rather than bulk processes i 2 - The 
primary energy conversion process in OPVs involves exci- 
tation of tightly bound electron-hole pairs (excitons) and 
charge carrier generation, i.e., a separation of charges 
across a heterointerface formed by an electron donor and 
acceptor (e.g., a conjugated polymer and a functionalized 
fullerene^). Because of the large exciton binding energies, 
the OPV system must provide a suitable energetic driv- 
ing force for the excitons to dissociate, making charge 
separation a key process that limits OPV efficiency. Ow- 
ing to their low cost and scalability, photovoltaics based 
on polymer and small molecule organics have received 
increasing attention in recent years 4^ However, basic 
understanding of the fundamental processes at the in- 
terfaces, such as charge separation and charge transfer, 
continues to lag behind that for inorganic PVs. 

An attractive alternative to all-organic solar cells is a 
hybrid solar cell that mixes an organic/polymeric sys- 
tem with an inorganic semiconductor/metal oxide such 
as CdSe^ Ti02r or ZnO^ as the electron acceptor. These 
inorganics are suitable for photovoltaic applications ow- 
ing to their relatively high electron affinity and mobility. 
The flexibility to pattern them into nanocrystalline struc- 
tures makes them desirable for highly ordered nanostruc- 
tured photovoltaic cells. 9 There have been many reports 
of hybrids combining poly(3hexylthiophene)(P3HT) and 
ZnO. 10 Their efficiencies, however, are much lower than 
those of fullerene based organic blends. Identifying the 
factors limiting the efficiencies of organic/inorganic solar 
cells remains a challenge. 

Owing to their highly ordered and controllable crys- 
talline structures, oligomers such as sexithiophene (6T) 
offer an ideal system to study the electronic structure 
at organic/inorganic interfaces^ - — Blumstengel et al. 
have recently reported atomic force and ultraviolet pho- 
toelectron spectroscopy (UPS) of 6T on ZnO surface 
showing that the heterostructure holds promise for pho- 
tovoltaic applications^ Depending on the surface metal- 



licity, directions, and coverage, the 6T molecules can take 
either a flat "lying-down" or an upright "end-on" ori- 
entation^ - — How the molecular orientations affect the 
interfacial electronic structures is not yet completely un- 
derstood. 

Here we carry out a first principles study of the sex- 
ithiophene and ZnO interface using hybrid functional 
density- functional theory (DFT). Our calculations ad- 
dress the question of orientation dependence and show 
that both the upright and fiat-lying orientations give rise 
to type-II interfacial band alignment that is necessary for 
charge separation. Their interfacial energy offsets are, 
however, very different. Controlling the orientations of 
organic oligomers on inorganic substrates offers an op- 
portunity for molecular engineering of the interface so as 
to achieve good photovoltaic performance. 

We carry out calculations using a plane-wave pseu- 
dopotential approach to DFT implemented in the VASP 
codei£ with PAW potentials H To study electronic prop- 
erties we employ the PBEO hybrid- functiona l 18 ! 19 in 
which a fraction of the nonlocal exact exchange energy 
is included in the exchange-correlation potentials. Al- 
though not completely free of the self-interaction error ^ 
hybrid functionals have yielded much improved band 
gaps for semiconductors including ZnO*2i - — Computa- 
tional methods that minimize the error in bandgap cal- 
culations are essential for accurate theoretical predictions 
of interfacial energy alignment involving organic or inor- 
ganic systems. 24 26 We use 6x6x6, T-point, and 1x2x1 
for fc-point sampling of bulk ZnO and the two super cells, 
respectively, and a plane-wave energy cutoff of 400 eV. 

Using the PBE generalized gradient approximation, 27 
we obtain lattice constants of a = 3.29A, c = 5.30A (exp. 
a = 3.25A, c = 5.21A 28 ) that agree with previous stud- 
ies 29 and a band gap of 0.7 eV for bulk ZnO. The PBEO 
hybrid functional gives rise to band gaps of 3.08 eV and 
3.16 eV for the bulk and the (1010) surface, respectively, 
in much better agreement with the experimental 3.37 
eV*2£ For the 6T molecule, we have obtained an energy 
gap of 1.57 eV in PBE and 2.7 eV in PBEO calculations 
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FIG. 1: (color online) Atomic structure of the 6T/ZnO (1010) 
interface where the 6T molecule is in (a) a flat-lying (N) and 
(b) an upright standing (L) orientation, (c) The interfacial 
structure after the relaxation for the N interface. The dashed 
arrows mark the average distance of the molecule above the 
ZnO trench rows. 



(at the PBE geometry), the latter compares well with 
the experimental value 30 of 2.8 eV and an earlier B3LYP 
hybrid functional calculation^ 

Although the growth of 6T on both nonpolar and polar 
ZnO surfaces has been demonstrated, 15 the structure of 
polar ZnO surfaces, and in particular the charge passiva- 
tion and surface reconstruction mechanism, are still un- 
der debate i 32 i 33 Here we consider the uncharged nonpolar 
(1010) surface that has the lowest cleavage energy and is 
easy to work with experimentally. 34 We have considered 
a ZnO (1010) surface slab including four ZnO double lay- 
ers (see Fig. [1]). Each surface unit cell contains a ZnO 
dimer row separated by a trench row along the [0001] 
direction. Surface relaxation has led to a strong tilt of 
the ZnO dimer as a result of the anions moving closer 
to the bulk, similar to what has been reported earlier^ 
The upright standing (L) and fiat-lying (N) interfaces 
are constructed by aligning the long (backbone) and nor- 
mal axis of the molecule, denoted by L and N, to the 
surface normal, as illustrated in Fig. [1] (a) and (b). Ow- 
ing to the incommensurate surface lattice parameters for 
6T crystal (a = 6.0 A and b = 7.9 A)^ 5 - and ZnO (1010) 
surface (a = 3.288 A and b = 5.303 A) along with a peri- 
odic boundary conditions, modeling the interface with a 
6T and matching ZnO slab would introduce an unrealis- 
tically large strain. Instead we have adopted a supercell 
that contains 2x2 and 10 x 2 ZnO surface unit cells 
for the L and N orientation, respectively. This neglects 
the intermolecular interactions. Each supercell contains 
one 6T molecule initially placed above the trench row. 
Recent calculations of Dag & Wang within the local den- 
sity approximation (LDA) have shown this position to be 
energetically favorable for a fiat-lying P3HT molecule on 
ZnO (1010) surface. 36 We include a vacuum of more than 
8 A along the interface direction to ensure isolation of 
the interface. 

Starting with a fully relaxed ZnO (1010) surface, we 



initially position the 6T molecule at a height of 2-3A (de- 
pending on the orientation) above the trench and relaxed 
the atomic positions of the molecule and the first two 
ZnO surface double layers while fixing position of the bot- 
tom ZnO layers. Upon the relaxation, the 6T molecule 
undergoes about 20° tilting relative to the surface in the 
N orientation due presumably to interaction between the 
molecule tt face and the tilted Zn-0 dimer s, as shown in 
Fig. [T](c), while remaining perpendicular to the surface 
in the L orientation. The final position of the molecule 
has an average adsorption distance of 2.8 and 1.7 A for 
the N and L orientations, respectively, above the surface. 
The former value is comparable to that of P3HT/Zn0. 36 
Because of the computational cost of the hybrid func- 
tional calculation (which is up to 30-50 times more ex- 
pensive than with the standard functionals) and a large 
supercell size, we have kept the interfacial geometries at 
the PBE level which has been shown to yield reasonable 
surface geometries^ Van der Waals (vdW) interactions 
which are absent in the PBE functional may influence 
the molecular tilting and adsorption distanced Tests in- 
dicate, however, that removing the tilting and reducing 
the adsorption distance by up to 0.3 A produce essen- 
tially the same electronic structure in our results. 

The adsorption energy of 6T on the ZnO interface is 
calculated from the difference between the total energy 
of the interface and that of the fully relaxed ZnO surface 
and 6T molecule. Within the PBE functional, we have 
obtained an adsorption energy of ~ leV for the N and 0.1 
eV for the L orientation. (The PBE0 functional gives 2 
eV and 0.4 eV, respectively, at the PBE geometry.) The 
former value seems much smaller than the binding energy 
of 2.8-2.9 eV in ref. [H. Owing to the lack of vdW inter- 
actions, neither PBE (PBE0) nor LDA functional gives 
an accurate account of the adsorption binding energy for 
physisorbed molecular interfaces^ An explicit inclusion 
of the van der Waals (vdW) interaction is, however, com- 
putationally demanding and well tested parameters are 
lacking for systems involving metal oxides. 38,39 Here we 
focus on the relative stabilities of different molecular ori- 
entations. Compared to the upright standing orienta- 
tion, the flat-lying one is much more energetically stable 
owing to a stronger orbital overlap between the molec- 
ular 7r wavefunctions and the substrate surface. This 
relative stability applies to the situation in which the 
molecule-substrate interaction is dominant (e.g., in a low 
coverage interface), but but the upright standing orien- 
tation would be stabilized if intermolecular interactions 
are taken into account. 

We apply the PBE0 hybrid functional to study the 
interfacial electronic structures. In Fig. [2fa) we plot 
the projected density of states (PDOS) of the upright 
standing interface. The Zn 3d band lies 6 eV below, 
and is well separated from the top of the valence bands 
which are composed of hybridized Zn 3d and O 2p states. 
The bottom of the conduction band (Zn 4s) lies 3 eV 
above the valence band maximum (VBM). The 6T states 
are discrete, with the highest occupied molecular orbital 
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FIG. 2: (color online) Projected density of states on the ZnO 
(blue solid lines) and 6T orbitals (red solid lines shaded from 
below) of the upright standing interface obtained from the 
(a) PBEO hybrid and (b) PBE functional. The dashed line 
denotes the Fermi energy which is at the same level of the 6T 
HOMO. 



(HOMO) level lying in the ZnO gap and the lowest un- 
occupied molecular orbital (LUMO) level about 2.5 eV 
above the Fermi level overlapping with the ZnO conduc- 
tion bands. Fig. EJb) plots the PDOS calculated from 
the PBE functional. The Zn 3d states are much shal- 
lower in energy, this pushes up the O 2p states and leads 
to a much smaller gap in ZnO compared to the case of 
PBEO. The ZnO conduction band minimum (CBM) also 
lies much lower in energy than in the PBEO case and over- 
laps with the 6T HOMO at the Fermi level. The band 
offset predicted using the PBE functional is unphysical. 
Henceforth we focus on the PBEO results. 

A suitable energy offset at the donor and acceptor in- 
terface is a key aspect for charge separation and charge 
transfer in OP Vs. To calculate the level alignment at 
the 6T/ZnO interface, we focus on the PDOS near the 
HOMO and LUMO. For the upright standing (L) inter- 
face as shown in Fig. 02 (a), the 6T HOMO and LUMO 
levels are both higher in energy than the ZnO VBM and 
CBM respectively, indicating that the interfacial align- 
ment is a staggered type-II interface ideal for organic 
photovoltaics. By comparing the energy positions, we 
find that the A^homo/vbm and A£ LUMO /cbm offsets 
are of 2.5 eV and 2 eV. From the UPS measurement, 15 
Blumstengel et al. have reported a A£? HO mo/vbm of 2.15 
eV and a A£ LUMO /cbm of 1.5 eV for 6T/ZnO (0001). 
The latter offset was obtained by adding the optical en- 
ergy gap, which take into account of the exciton binding 
energy of 6T and ZnO, to the occupied levels. Using the 
same optical gaps, we have obtained a A£? L umo/cbm off- 
set of 1.8 eV. Despite the absence of intermolecular inter- 
actions, our data compares well with the UPS result and 
confirms that the experimental band offsets arise from 
upright standing molecules. The small differences maybe 
attributed to the surface dipole on ZnO (0001) which af- 
fects the band offset. Compared to the exciton binding 
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FIG. 3: (color online) Projected density of states for the (a) 
L oriented and (b) N oriented interface. The vertical dashed 
line marks the 6T HOMO and LUMO and the ZnO VBM 
and CBM. The LUMO and CBM are obtained by adding the 
corresponding optical gap for 6T (2.7 eV) and ZnO (3.35 eV) 
to the occupied levels. Below the panels are calculated energy 
offsets for the 6T/ZnO interfaces. 



energy of 6T molecules (~ 0.4 eV), the LUMO/CBM off- 
set is significantly higher, indicating a sufficient driving 
force for exciton separation across this interface. 

The flat-lying (TV) interface, as shown in Fig.[3fb), also 
offers a type-II interface. However, the HOMO/VBM 
and LUMO/CBM offsets, 1.4 eV and 0.75 eV (0.84 eV 
from PDOS), respectively, are almost 50% lower than 
those for the L interface, albeit still sufficient for charge 
separation. The percentage change of the offsets with 
respect to the molecular orientaion compares very well 
in magnitude with the recent experimentally measured 
values for sexiphenyl (6P)/ZnO interfaces^ We further 
extract the open circuit voltage (V oc ) that plays a sig- 
niciant role in photovoltaic efficiency. This is defined as 
the difference between the quasi-Fermi levels of the sep- 
arated holes and electrons and can be approximated by 
the offset between the HOMO of the donor and the CBM 
of the acceptor. For the N oriented interface, we find a 
V oc of 1.9 eV which is significantly higher than 0.9 eV for 
the L oriented interface. Our results suggest a sensitive 
molecular orientation dependence of the energy offsets 
at hybrid interfaces. A control of the molecular orien- 
tations should provide an effective approach to increase 
the output potential and the OPV efficiency. 

To understand the orientation dependence of the en- 
ergy offsets, we carried out an analysis of the frontier 
orbitals of the 6T molecule at the upright (L) and fiat- 
lying (N) interfaces. For both cases, the 6T HOMO is 
localized in the ZnO band gap while the LUMO overlaps 
with the ZnO conduction band in energy. This over- 
lap has important implications in the electron transfer 
across the interface. Fig. H(a) shows the spatial distri- 
bution of the electron density of the LUMO at the up- 
right interface. The LUMO is a 7r* orbital and is almost 
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FIG. 4: (olor online) Isosurface of the electron density (the 
isovalue equals to the average of density) of the 6T LUMO 
on the (a) upright and (b) flat-lying interface. Panel (a) only 
shows half the number of thiophene rings in 6T. 

completely localized on the molecule chain and displays 
little electronic coupling with the ZnO substrate despite 
the overlap in energy. In contrast, the LUMO orbital at 
the flat-lying interface, as shown in Fig. HJb), clelocalizes 
over both the molecule and the substrate (especially on 
the molecular side closer to the surface because of the tilt 
of the thiophene rings against the surface), indicating a 
strong electronic coupling between the molecule and the 
substrate at the interface. The fact that the electronic 
coupling is much stronger for the flat-lying case com- 
pared to the upright case is consistent with the higher 
adsorption energy we find for the N orientation. We 
can also understand the orientation dependence of the 
interfacial energy offsets and value of V oc qualitatively 
from the strength of the interfacial electronic coupling. 
In case of the flat-lying interface, the molecular orbital 



levels have shifted significantly closer to the substrate 
levels as a result of the stronger coupling. This signifi- 
cantly reduces the LUMO/CBM and HOMO/VBM offset 
between the two materials while raising the V oc . 

In summary, we have employed a hybrid density 
functional to study the electronic structure of sexithio- 
phene/ZnO interfaces. The upright and flat-lying ori- 
entations both offer suitable type-II interface that are 
ideal for charge separation across the interface. However 
the magnitudes of the energy offsets at the interface are 
strongly dependent on the molecular orientations. The 
calculated offsets for the upright orientation agree well 
with the experimental values obtained from UPS^ and 
the relative change of the offsets with respect to the ori- 
entations also compares fairly well with the measured 
change for 6P/ZnO interfaces^ Our analysis suggests 
that the varied strength of the electronic coupling be- 
tween the molecule and the substrate should account for 
the orientation dependence of the energy offsets. 
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